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ABSTRACT 

We present a morphology study of intermediate-redshift (0.2< z <1.2) luminous infrared galaxies 
(LIRGs) and general field galaxies in the GOODS fields using a revised asymmetry measurement 
method optimized for deep fields. By taking careful account of the importance of the underlying 
sky-background structures, our new method does not suffer from systematic bias and offers small un- 
certainties. By redshifting local LIRGs and low-redshift GOODS galaxies to different higher redshifts, 
we have found that the redshift dependence of the galaxy asymmetry due to surface-brightness dim- 
ming is a function of the asymmetry itself, with larger corrections for more asymmetric objects. By 
applying redshift-, IR-luminosity- and optical-brightness-dependent asymmetry corrections, we have 
found that intermediate-redshift LIRGs generally show highly asymmetric morphologies, with implied 
merger fractions ~ 50% up to z=1.2, although they are slightly more symmetric than local LIRGs. 
For general field galaxies, we find an almost constant relatively high merger fraction (20-30%). The 
£>-band LFs of galaxy mergers are derived at different redshifts up to z=1.2 and confirm the weak 
evolution of the merger fraction after breaking the luminosity-density degeneracy. The IR luminosity 
functions (LFs) of galaxy mergers are also derived, indicating a larger merger fraction at higher IR 
luminosity The integral of the merger IR LFs indicates a dramatic evolution of the merger-induced IR 
energy density [(l+z)~( 5 ~ 6 )], and that galaxy mergers start to dominate the cosmic IR energy density 
at z > ~1. 

Subject headings: infrared: galaxies - galaxies: interactions 



1. INTRODUCTION 

Our understanding of galaxy formation and evolution 
has advanced dramatically in the past decade, with well 
determined cosmic evolution of the comoving star for- 
mation density, g alaxy stellar mass and galaxy metal- 
licity content (see lCowie fc Bargerll2008l and references 
therein). What mechanism drives this evolution? Galaxy 
interactions and mergers play a key role in the theory of 
galaxy evolution, through transforming the galaxy mor- 
phologies, inducing violent st arbursts, and feeding the 
central massive black holes (jMihos fc Hernquistl fl996: 
iSpringel et al.ll2~005f) . However, it is still unclear from 
observations to what extent galaxy mergers actually play 
the roles predicted for them. 

Galaxy mergers can be found from observations ei- 
ther morphologically or kinematically. Morphologically- 
identified mergers include galaxies with tidal tails, wisps 
and/or multiple nuclei. They can be found through ei- 
ther visual classifications or quantitative measurements. 
While visual classification is the testbed to develop the 
quantitative methods, it is time-consuming to classify 
the enormous numbers of galaxies found in deep surveys. 
Moreover, distant galaxies can change their apparent 
morphologies due to surface-brightness (SB) dimming. 
Such an effect can be corrected better through quanti- 
tative measurements. Several quantitative morphology 
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techn i ques have been developed, CAS (|Abraham et al.l 
11994 119961: IConselice et all I2000L 12003ft and Gini-M 20 
(| Abraham et al.ll2003l : lLotz et al.ll2004ft . They classify 
mergers by identifying galaxies in a pre-defined parame- 
ter space. Although the interpretation of the parameter 
space is somewhat uncertain, studies based on the same 
quantitative definition can be compared and are less sub- 
jective than visual classification. Galaxy mergers can 
also be identified through kinematical pairs of galaxies 
or galaxies with complicated internal velocity fields (i.e. 
neither pressure nor rotationally supported). Identify- 
ing them requires time-consuming spectroscopic observa- 
tions, which h ave been conducted for only a limited num - 
ber of objects (|de Ravel et al.ll2008t TNeichel et al.ll2008h . 

Different approaches to merger identification are sen- 
sitive to different stages of the galaxy merging process. 
Pair-identification algorithms find separated interacting 
galaxy pairs while the morphologically-based algorithms 
usually identify galaxies during the first pass and final 
coalescence where the galaxy morph ology is significantl y 
disturbed by gravitational torques (jLotz et al.l I2008b1) . 
Galaxy mergers with complicated interval velocity fields 
seem to span a long er timescale of the merging process 
(|Neichel et al.l [20081 ) . Nevertheless, no method can iden- 
tify all mergers. In addition, a fraction of galaxies identi- 
fied as mergers by any method are not necessarily major 
mergers. For example, minor mergers of gas-rich galaxies 
can also have highly disturbed morphologies (Lotz et al. 
2009 in preparation). 

Morphology studies of high-redshift galaxies show that 
the Hubble sequence is in place by z ~1 and that high- 
redshift galaxies ar e associated more freque ntly with pe- 
culiar features (e.g. Brinch mann et al.lll998l ). Some stud- 
ies have found that the galaxy merger fraction shows 
strong redshift evolution, characterized by (1 + z) m with 
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m > 2 (e.g. iLe Fevre et all I200C 


; iPatton et all 120021: 


Conselice et al.ll2003l ICassata et al. 


120051: Ide Ravel et al.l 


2008}) , implying that the cosmic 


star formation his- 



tory (SFH) is at least partly driven by galaxy merg- 
ers. Other studies, however, found a slow evolution 
(e.g JBundv et all l2004t iLin et al.l 120081: Ide Ravel et all 
I2008D or even no evolution (|Lotz et al.ll2008a|) . The 
discrepancies among different studies may be caused 
by l ow-number statistics , morphological K-corrcctions 
(e.g. iPapovich et al.ll2003l) . the method to identify merg- 
ers (|Lotz et al.l l2008bT ) and the physi cal properties (e.g . 
gas fractions) of merger progenitors (jLotz et al.l l2008bL 
Nevertheless, almost all studies have found relatively low 
merger fractions (<20%) at z < 1. 

A more direct way to constrain the role of galaxy 
mergers in the cosmic SFH is to investigate morpholo- 
gies of luminous infrared galaxies (LIRGs, Lir > 10 11 
L Q ). LIRGs show strong redshift evolution and start 
to dominate the cosm i c IR energy density at z > 0.7 
(|Le Floc'h et al.l 120051 : iPerez-Gonzalez et al.l l2005| ). In 
the local universe, LIRGs are mainly (~50% ) triggered 
by major mergers ( jSanders fc Mirabel 1 1996tl . However, 
morphology studies of intermediate-re dshift LIRGs indi- 
cate low merger f ractions (~10-30%) (JZheng et al.ll2004| : 
Bell et al.l 120051: iBridge et all 120071 : fLotz et al l l2008a[ 
Melb ourne et al I I2008D . These findings have led to 
claims that the cosmic SFH is driven by some less vi- 
olent mechanisms, such as accretion or gas consumption 
(jNoeske et al.ll2007t iDaddi et al.ll200l . 

In IShi et all (|2006f ) , we measured the galaxy asymme- 
try for LIRGs in the Ultra Deep Field (UDF), which pro- 
vides limiting SB for galaxies at z=\ comparable to that 
often obtained for local ones (hb = ~25.3 mag arcsec -2 
at lOcr in the AB magnitude system). The merger frac- 
tion obtained in this study is several times higher than 
others, 40±24% and 26±10 % for LIRGs and general 
field galaxies (Mb < -19.25) at z ^0.7, respectively. The 
comparison (see Figure QJ between asymmetry measure- 
ments using the GOODS images and ones using the UDF 
images for the same galaxies indicates substantial deficits 
in the indicated asymmetry in the shallower GOODS 
images. These two results motivate us to suspect that 
other studies based on relatively shallow imaging (most 
of which are shallower than GOODS) may underestimate 
significantly the role of SB dimming in galaxy morphol- 
ogy measurements. If so, the role of mergers in triggering 
luminous episodes of star formation may be seriously un- 
derestimated. 

To test the result of lShi et all (|2006l ) with much higher 
statistical significance and better constrain the role of 
galaxy mergers in the cosmic SFH, in this paper we carry 
out detailed asymmetry measurements and corrections 
for all galaxies (~16,000) with m z < 25 including ~7,500 
galaxies and ~1000 LIRGs at z < 1.2 in the GOODS 
field. In § [21 we describe our revised asymmetry esti- 
mate method. In § [3l data for a complete local LIRG 
sample and the GOODS ACS-, redshift- and Spitzer- 
observations are presented, as well as the concentration 
and asymmetry measurements for both local LIRGs and 
GOODS galaxies. §|4]shows the evolution of the observed 
asymmetries of GOODS galaxies (§ I4.1j) and the results 
after asymmetry corrections, including the redshift evo- 
lution of the merger fraction in LIRGs (§ |4.2|1 . the evolu- 
tion of the merger fraction in general field galaxies (§ 14. 3| 
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Fig. 1. — Comparison of galaxy morphologies using the GOODS 
z-band images and the UDF z-band images for about 250 UDF 
objects with relatively high S/N and large size. The median er- 
ror bar is shown in bottom-right corner. The solid line is for 
Asy(GOODS)=Asy(UDF). 

and the infrared and rest-frame B-band luminosity func- 
tions of galaxy mergers (§ I4.4[) . In § O we discuss the 
role of mergers in galaxy star formation activity since 
z ~1. The conclusions are in § [6] The technical details 
are given in § [7] for our revised asymmetry method and 
m § HI f° r asymmetry deficits of redshifted local LIRGs 
and low-z (z=0.2-0.4) GOODS galaxies due to SB dim- 
ming. Throughout the paper, we adopt a cosmology with 
iJ =70 km s^ 1 Mpc" 1 , ft m =0.3 and O\=0.7. All mag- 
nitudes are defined in the AB magnitude system. 

2. GALAXY STRUCTURE MEASUREMENT 

2.1. Galaxy Size And Concentration 

Quantitative galaxy morphologies are always measured 
within well-defined apertures in order to probe the same 
physical regions of galaxies; these apertures are opti- 
mized to include as much galaxy flux as possible while 
minimizing the effect of noise plus SB dimming. When 
comparing low- to high-rcdshift galaxy morphologies, SB 
dimming (SB oc (1 + z) 4 ) can artificially decrease the 
galaxy size. To avoid such an effect, the galaxy size can 
be measured through the dimcnsionless parameter: 



V (r) = /(r)/(J(r)> 



(1) 



where I(r) is the surface brightness at radius r and (I(r)) 
is the mean surface brightness wit hin radius r ([Petrosianl 
1976 ). In the standard approach (|Conselice et al.ll2000l 
20031) . the Petrosian radius R p is defined at ?7(r)=0.2 
and the galaxy aperture radius is defined to be 1.5i? p . 
The surface brightness I(r) can be measured within ei- 
ther circular apertures or elliptical apertures. For in- 
clined non-interacting galaxies, the mean galaxy elliptic- 
ity and position angle trace the true galaxy light well. 
The R p defined within an elliptical aperture is gener- 
ally larger than that defined within a circular aperture 
and can be two times larger for extremely inclined sys- 
tems. However, we have found that interacting galaxies 
usually have large variations in their position angles and 
ellipticities. The measured galaxy size for them depends 
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on the adopted position angle and ellipticity. Because 
of this uncertainty for the interacting systems, we adopt 
circular apertures universally for all galaxies in this pa- 
per. For most local LIRG systems, we found that the 
difference in i? p between circular apertures and ellipti- 
cal apertures with mean ellipticity and position angle is 
within 20%. Therefore, using circular apertures does not 
introduce large errors. 

The concentration parameter measures how compact 
the light distribution is. It has been shown to correlate 
well with the Hubble type (jKentl 119851 : [Bershad v et al.l 
2000). Early-type galaxies generally have more compact 
light distributions than do late- type ones. With the curve 
of growth measured within the galaxy aperture, the con- 
centration can be defined as the ratio of the two radii 
enclos ing t wo fixed fraction s of th e light. Following iKentl 
(119851 ) and lBershadv et alj (|2000f ). we have used the fol- 
lowing definition in this paper: 
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where i?so and R20 are the radii that enclose 80% and 
20% of the total light, respectively. 



2.2. Revised Asymmetry Parameter 
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param eter (^4) was first introduced 
and subsequently developed 
120031 ). It quantitatively de- 
scribes the level of the galaxy asymmetry and provides 
a direct measurement of the importance of asymmetric 
sub-structures in galaxy morphologies, such as merging 
companions and tidal tails, which are used to identify 
mergers in visual classification. Accurate measurements 
of the asymmetry parameter are thus critical to under- 
stand the role of galaxy mergers in galaxy evolution. The 
galaxy asymmetry is defined as: 



X/ I ^galaxy +e 



rl80° 
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where ^4 ga iaxy+noisc is the asymmetry of the galaxy signal 
plus the noise, igaiaxy+noisc is the image of galaxy signal 
plus the noise and I^i° Ky+noisc is /gaiaxy+noiso with the 
image rotated by 180° around a rotation center. The 
true pure galaxy asymmetry is then given by: 



^-galaxy fhl^^galaxy+noisc) A T 



(4) 



where min(A ga iaxy+noisc) is the minimum of the galaxy 
signal plus noise asymmetry, and A™ jg C is the noise cor- 
rection obtained by rotating a background image around 
a center and normalizing by the galaxy brightness: 



a corr 
noise 
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galaxy+noisc 



(5) 



where B is the background image without any object. 

The overall measurements of galaxy concentration and 
asymmetry are composed of the following steps: 

(1) Guess the initial center. 

(2) Measure the R p and concentration. 

(3) Search for the rotation center that gives the minimum 
of ^4gaiaxy+noise within a 1.5i? p aperture. 

(4) Use the new rotation center to re-measure the R p and 
concentration. 



(5) Use the above rotation center and new R p to search 
for a new rotation center that gives the minimum of 
^gaiaxy+noiso within a 1.5i? p aperture. 

(6) Correct for the noise ^noiL- 

The galaxy centers can be easily located for non- 
interacting galaxies but not for interacting ones. Here, 
we adopt the rotation center giving the minimum of 
^.gaiaxy+noise to measure the galaxy size as is done in 
steps 1-3. Note that the asymmetry rotation center 
in step 5 is usually not very different from that in 
step 3, as the "walk ing around" method invented by 
IConselice et al.l (|2000l ) is generally robust to find the min- 
imum A ga i axy+no i se . For discussions about the galaxy 
size definition, the algorithm searching for the minimum 
^galaxy +noise, and the dependence of - A ga i axv + noise on 
the resoluti on and correlated noise, see IConselice et al.1 
(l2000ll2003l) . 

The technical details about our new A£°g e are given 
in § Here we present a summary of the procedure 
to determine this parameter. A set of 1000 randomly- 
produced noise asymmetry measurements is carried out 
by putting circular regions in the background image 
around the target. The value of this distribution at the 
15% probability low-end tail is used as ^4 n °" c , the me- 
dian of true noise corrections. The error (~3cr) in the fi- 
nal measured galaxy asymmetry is taken to be two times 
the standard deviation of these randomly-produced noise 
asymmetries. In reality, some galaxies are always present 
in the field of targets. To account for this problem, a suc- 
cess rate is defined for the one-thousand circular region 
placements as the fraction of circular regions containing 
no galaxy signal indicated by the SExtractor segmenta- 
tion image. Circular regions containing any galaxy signal 
are not used. More sets of one-thousand placements are 
generated until one thousand successful measurements 
are reached. If the success rate for one set of 1000 place- 
ments is lower than 50%, the circular region size for the 
following set is taken to be 80% of this set. Then the 
measured noise asymmetry is rescaled to that with the 
original size by assuming that the noise asymmetry is 
proportional to the aperture area. 



3. GALAXY SAMPLE 

3.1. Local LIRG Sample 

Local LIRGs and ULIRGs are generally galaxies un- 
dergoing mergers with morphological signatures of tidal 
tails, multiple nuclei an d highly asymmetric features 
(|Sanders fc Mirabel! [i~996l ). To account for the redshift- 
dependence of galaxy merger morphologies due to SB 
dimming, we measured the morphologies of local LIRGs 
and of local LIRGs redshifted to different redshifts. We 
retrieved HST ACS /WFC-F435W images of 88 local 
LIRGs observed in Program ID - 10592 (PI - Aaron 
Evans). This local LIRG sample is a complete sub- 
sample of the IRAS Revise d Bright Galaxy S ample 
(RBGS: i.e., / 60Mm > 5.24 Jv: ISanders et al.ll2003f ) above 
Lik = 10 114 L . The sample size is large enough for sta- 
tistically valid comparisons. The distance of this sample 
covers the range from 35 to ^350 Mpc. The median dis- 
tance of 135 Mpc corresponds to a physical resolution of 
80 pc. These galaxies are also bright enough to be useful 
when redshifted to higher redshifts. These local LIRGs 
are currently experiencing merging processes, spanning a 
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wide range of merging stages from well-separated galaxy 
pairs to the final merging remnants. 

To obtain quantitative morphology measurements, 
foreground stars and background galaxies should be re- 
moved from the images first. Given the importance 
of the background region in accurate morphology mea- 
surements (see § E]), we carried out a careful removal 
of contaminators, especially given that some of the 
LIR Gs have tens of foregrou nd stars. We used SExtrac- 
tor (|Bertin fc Arnouti [l996) to obtain the segmentation 
map of each image. To make sure the extended low-SB 
emission was included in the object, we first rebinned the 
image by a factor of 4x4 and then adjusted the parame- 
ters in SExtractor so the extended emission was detected. 
This segmentation map was then resampled to the orig- 
inal resolution by assuming that the original pixels be- 
longing to a rebinned pixel conserve the segmentation 
value. We then cleaned the neighbourhood of the target 
by replacing pixels of non-target object pixels randomly 
with values of background pixels. As the neighbourhood 
is dominated by the background pixels (> 95%), such re- 
placement does not introduce any significant correlation 
in the final cleaned neighbourhood background. Some- 
times a target was identified as several different objects 
by SExtractor. We visually inspected each image to make 
sure such substructures were still included in the target 
during the neighbourhood clean. Sometimes foreground 
stars lie within the target region and cannot be separated 
by SExtractor. We used the IRAF package IMEDIT to 
remove these contaminators further. 

As the local LIRG sample is a flux-limited sam- 
ple instead of a volume limited sample, we have ap- 
plied weights of y^— /^2y^— to each object to obtain 
an equivalent volume-limited asymmetry distribution, 
where V max is measured at the redshift where the object 
with a given IR luminosity has feoum = 5.24 Jy. 

3.2. Galaxies in the GOODS Field 
3.2.1. HST Images and Morphology Measurements 

The GOODS field consists of two sub-fields, GOODS- 
North (GOODS-N) and GOODS-South (GOODS-S), im- 
aged with HST/ ACS in four filters B (F435W), V 
(F606W), i (F775W) and z (F850LP) (|Oiavalisco et all 
[200l M. Giavalisco and the GOODS Team, 2008, in 
preparation). The field centers (J2000.0) are 12h36m55s, 
62°14'15" for the GOODS-N and 3h32m30s, -27°48'20" 
for the GOODS-S. The survey area is 320 arcmin 2 
with BViz-b&nd coverage and 365 arcmin 2 with Viz- 
band coverage for each field. GOODS Version 2.0 pro- 
vides exposure times of 7200, 5650, 8530 and 24760 
sees in GOODS-N and 7200, 5450, 7028 and 18232 
sees in GOODS-S. The final image has a pixel scale of 
0.03"/pixel. The GOODS-S and GOODS-N fields are di- 
vided into 18 and 17 sections, respectively. An individual 
image file is released for each section. 

Object s were detected in z-ban d with the SExtractor 
package (jBertin fc Arnoutslll996h and photometry was 
measured in all four bands. We produced the segmen- 
tation images that define the galaxy pixels using the re- 
leased parameter files of SExtractor plus weight map im- 
ages, and final science images. We used galaxy magni- 
tudes defined by SExtractor MAG_AUTO. 

We carried out morphology measurements in the 



Viz-bands for 16,708 objects with m z < 25, SEx- 
tractor CLASS_STAR < 0.9, FLUX.RADIUS1 < 100, 
IMAFLAGSJSO < 16 and not within 33 pixels 
of the edge of each SECTION field. SExtractor 
FLUX_RADIUS1 is the radius in pixels enclosing 20% 
of the flux and FLUX_RADIUS1 < 100 excludes 9 artifi- 
cial objects that are long narrow bright belts across the 
field. IMAFLAGSJSO < 16 excludes objects within 33 
pixels of the field edge. The morphology measurement 
was first carried out in z-band starting with cutting out 
science and segmentation images for each object. As the 
galaxy radius is defined to be l.5R p and the asymmetry 
uncertainties are measured by putting circular regions 
randomly in background regions around the objects, the 
image size was cutout with a size of 9R p x9R p . Then the 
z-band concentration and asymmetry were measured as 
described in § 12.21 Briefly, the residual sky was first 
subtracted using the mean values of all pixels with zero 
segmentation values. Companions were defined as all 
pixels with non-zero segmentation values not equal to 
the target's value. They were removed by replacing their 
pixel values randomly with those of sky pixels. Note 
that the removed companions are those well separated 
from targets and thus do not contribute to the target 
asymmetries. The initial center was given by the astro- 
metrical position of the object. The first set of concen- 
tration/asymmetry measurements gave new centers and 
the second measurements were carried out to give R p , 
concentration, and asymmetry values. The asymmetry 
uncertainty was measured in 1000 circles randomly put 
in the sky region, excluding object and companions (for 
details, see the last paragraph of § [Tj) . To evaluate the de- 
pendence of the asymmetry uncertainties on the S/N, two 
types of S/N are defined. The total S/N is the total sig- 
nal divided by the total sky noise for galaxy pixels within 
the l.5R p radii. Not e that sky and com panion pixels are 
excluded. Similar to lLotz et all (|2008af ). (S/N) is defined 
to be the arithmetic (not quadratic) average of the S/N of 
each galaxy pixel. The asymmetry/concentration mea- 
surements failed for 212 objects (175 of them are stars; 
the remaining 37 objects either have extremely low SB 
or are near the edge, within 50 pixels). No systematic 
correction is made for these objects, as the SB limiting 
cut will be applied in the final morphology catalog and 
they are a very small fraction (0.2%) of the total number 
of objects. The I^i-band concentration and asymmetry 
were then measured for all objects with successful z-band 
concentration and asymmetry measurements. 

3.2.2. Spitzer MIPS and IRAC Data 

MIPS and IRAC data for the GOODS-N and GOODS- 
S fields were drawn from the Spitzer GTO observa- 
tions of the larger area covering each field, 1.5° xO. 5° 
for MIPS and 1.0°x0.5° for IRAC. The reduction of 
the 24 /xm ima ges was carried out with the MIPS Data 
Analysis Tool ([Gordon et~aT1l2005f l The detailed data 
reduction, object detection and photometr y mea sure- 
ment procedures are given i n iPapovich et all (|2004l ) and 
iPerez-Gonzalez et alJ (|2005[ ). The final MIPS 24 ^m cat- 
alog is 50% complete at 60 /zJy. 

3.2.3. Redshifts And Morphology- Redshift Catalog 

For the GOODS-S field, the spectroscopic redshifts 
were obtained from Version 3.0 of the FORS2 catalog 
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Fig. 2. — The S/N vs. asymmetry uncertainty (3<r), where the 
upper panel shows the mean S/N and the lower panel is for the total 
S/N. The S/N is only correlated with the asymmetry uncertainty 
for galaxies with S/N > 30 and (S/N) > 3. 

(IVanzella et alJl2008D and Version 1.0 of the VIMOS cat- 
alog (jPopesso et al.lr2008f ). Only solid redshifts with red- 
shift quality z q —A! were used. For GOODS-N, we used 
spectroscopic redshifts with z q —3 or 4 from the Team 
Keck Treasury Redshift Survey (TKRS) catalog. For 
objects with out secure spectroscopic red shifts, we used 
the catalog of lPerez-Gonzalez et al.l (|2008f ) . who obtained 
photo-z based on photometry covering from UV to IRAC 
bands. The photo-z error a z /(l+z) is < 0.2 for 95% of 
the redshifts and < 0.1 for 88% of the redshifts. The me- 
dian a z /(l+z) is 0.03. The photo-z error is small enough 
for our purpose, i.e., to determine the rest-frame band of 
the galaxy image. These redshift catalogs were matched 
to the GOODS z-band morphology catalog using a search 
radius of 0.5". The redshift was assigned to the nearest 
one if multiple objects were present within a search ra- 
dius. The redshift completeness for m z < 25 is 68%; it 
is mostly determined by the IRA C detections for which 
the ph oto-z can be calculated in lPerez-G onzalc z~et al.l 
( 2008). We refer to the sample of objects with both mor- 
phology and redshift measurements as the morphology- 
redshift catalog. 

To account for the morphological K-correction, the red- 
shift evolution of the galaxy morphology was determined 
in the rest-frame B-band for objects at z < 1.2. The red- 
shift range of z < 1.2 was divided into five redshift bins, 
[0.2, 0.4], [0.4, 0.6], [0.6, 0.8], [0.8, 1.0] and [1.0, 1.2]. The 
rest-frame B-band morphologies in a given redshift bin 
were defined as the morphologies in the band nearest to 
the redshifted rest-frame B-band, i.e, the observed-frame 
V— , V—, i—, i— and z-band for the redshift bins from 
low to high. 

To obtain the optical counterparts of the MIPS sources, 
we first obtained the IRAC 8/xm counterparts using a 
search radius of 1". Then for MIPS sources with IRAC 
counterparts, the optical counterparts were defined as ob- 
jects in the morphology catalog within 1" of the IRAC 
counterparts. For MIPS sources without IRAC coun- 
terparts, optical counterparts were searched for within 
2.5" of the MIPS source. For multiple objects within 
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Fig. 3. — The galaxy radius vs. magnitude in u~, i-, z- bands 
for galaxies selected with m z < 25. Areas with different grayscales 
correspond to different completeness cuts of galaxies with reliable 
morphology and redshift measurements (see ^? 13.36 . The solid lines 
enclose the 70% complete area where the apparent surface bright- 
ness, size and magnitude are labelled. 

a search radius, the nearest one was defined to be the 
optical counterpart. For MIPS sources with known op- 
tical redshifts, the total 8-1000 /im IR luminosity was 
obtained based on the observed 24 flux density and 
the star-forming templates from lRieke et al.l (|2009f) . 

3.3. Reliable Morphology And Completeness Cut 

The redshift-morphology catalog was further limited 
to objects with reliable asymmetry measurements, de- 
fined as those with i? ga i=1.5i? p > 15 pixels (0.03"/pixel) 
and error(^L) < 0.1 in a given band. The accuracy of 
the galaxy asymmetry depends on the resolution and 
S/N. For a galaxy with small size, the resolution is not 
high enough to resolve galaxy str uctures and the asym - 
metry becomes artificially small ([Conselice et al.l [2000). 
For galaxies with low S/Ns, the asymmetry uncertainty 
is large due to the noise. For the study in this pa- 
per, we used an asymmetry uncertainty cut (error(A) 
< 0.1 at approximately 3cr level) instead of a S/N cut, 
as our revised asymmetry method can characterize the 
uncertainty due to complicated background structures. 
Figure [2] shows the advantage of the asymmetry uncer- 
tainty cut over the S/N cut. While both S/N and (S/N) 
are tightly correlated with the asymmetry uncertainty 
at high S/N (e.g., (S/N) > 3), the scatter of the cor- 
relation becomes large at (S/N) < 3. A high S/N cut 
(e.g., (S/N) > 3) will exclude many objects with rela- 
tively small asymmetry uncertainty (error (A) < 0.1) and 
a low S/N cut (e.g., (S/N) > 1) will introduce objects 
with large asymmetry uncertainty (error(A) > 0.1). Our 
error(A) < 0.1 criterion only excludes objects with large 
asymmetry uncertainties. Most objects with error(j4) > 
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Fig. 4. — The first three panels show the rest-frame B-band mag- 
nitude, the galaxy radius and the rest-frame B-band surface bright- 
ness as functions of redshift for all galaxies with m z < 25 (dots) 
and MlPS-detected galaxies (open circles). The heavy solid lines 
correspond to the 70% completeness cut of the redshift-morphology 
sample (see § I3,3j l, The last panel shows the total IR luminosity 
of MlPS-detected galaxies as a function of redshift where the solid 
lines correspond to the 50% completeness cut of the 24^tm detec- 
tions. 

0.1 have relatively inhomogeneous backgrounds. All ob- 
jects with error(A) > 0.3 are at the edge of the field where 
the exposure is relatively shallow and shows a large gra- 
dient toward the edge. 

We now have a redshift-morphology catalog with red- 
shift information and reliable morphology measurements. 
The redshift evolution of galaxy morphologies will be 
characterized in a subsample of this catalog that is com- 
plete to a certain value. Figure [3] shows the completeness 
cut at 50%, 70% and 90% of the redshift-morphology 
catalog in i? ga i vs. magnitude in three bands, where 
the completeness is defined as the ratio of the num- 
ber of galaxies with redshift measurements and secure 
morphologies to the total number of galaxies in a i? ga r 
magnitude bin. The absolute values at different redshift 
bins are given in Table [TJ Here, we do not account for 
incompleteness of the HST/ ACS catalog, as its incom- 
pleteness at m z < 25 is negligible. 

Figure 2] shows the absolute magnitude, galaxy size, 
surface brightness and IR luminosity as a function of 
redshift for all m z < 25 and MlPS-detected objects with 
redshift measurements. The absolute rest-frame _B-band 
magnitude is m easured through the version 4.1.4 KCOR- 
RECT code (|Blanton et all l200l iBlanton fc Roweisl 
[2001 using four ACS bands. For the first three quanti- 
ties, the 70%-complete redshift-morphology cut is drawn 
in each redshift bin and the 50%-complete IR flux limit 
is drawn for the IR luminosity. As shown in the figure, 
most objects with m z < 25 and z < 1.2 are included in 
the final redshift-morphology catalog. This is because 



the incompleteness of the redshift-morphology catalog is 
mainly caused by the redshift incompleteness instead of 
the criterion of secure morphologies. 

4. RESULT 

4.1. Observed Asymmetry Distributions of GOODS 
Galaxies 

Two main results described in the Appendix are im- 
portant for the study of morphology evolution. First, as 
shown in § [71 we found that the underlying structure of 
the sky background is important in accurate asymmetry 
measurements and we propose a revised asymmetry mea- 
surement approach. Secondly, by redshifting local LIRGs 
and low-z (z=0.2-0.4) GOODS galaxies to higher red- 
shifts, asymmetry deficits due to SB dimming for these 
two galaxy populations are derived as a function of red- 
shift. As shown in § 18.31 originally more asymmetric ob- 
jects show larger asymmetry deficits. Such dependence 
indicates that to recover the intrinsic distribution of a 
galaxy population at a given redshift, a low-z galaxy 
population with intrinsically similar asymmetry should 
be used to construct the asymmetry corrections. Such a 
low-z galaxy population can be defined through compar- 
ing the observed asymmetries of redshifted low-z galaxies 
to those of the galaxy population at a given redshift, as 
the observed asymmetries still on average correlate with 
original asymmetries even given the larger asymmetry 
deficits for originally more asymmetric objects (see Fig- 
ure [J). 

Figure [5] shows the observed asymmetry of local 
LIRGs, GOODS LIRGs and GOODS field galaxies at 
different redshifts compared to the observed asymmetry 
of redshifted local LIRGs and redshifted low-z (z=0.2- 
0.4) GOODS galaxies. The GOODS LIRGS are the ob- 
jects at iin > 2.5 xlO 11 Lq and satisfying the 70% com- 
pleteness cut of the redshift-morphology catalog at the 
corresponding redshifts (see Figure [3] and Table [1} . The 
GOODS field galaxies within a redshift interval are all 
galaxies satisfying the 70% completeness cut of the high- 
est redshift bin. For the redshifted local LIRGs and low-z 
GOODS galaxies, the observed asymmetry becomes pro- 
gressively smaller at higher redshift as caused by SB dim- 
ming (see Appendix for discussion). 

Figure [5] shows a trend of the observed asymmetry 
for different subsets of the galaxy population within 
a given redshift bin: ^(redshifted z=0.2-0.4 GOODS 
galaxies) » A(GOODS galaxies) < A(GOODS LIRGs) 
< A(redshifted local LIRGs). This observed asymmetry 
trend and Figure [1] (i.e., that the observed asymmetry on 
average correlates with intrinsic asymmetry and that the 
asymmetry correction also increases with intrinsic asym- 
metry), indicate that the asymmetry corrections based 
on the asymmetry deficits of redshifted low-z GOODS 
galaxies should give a reasonable estimate of the true 
corrections for GOODS field galaxies. The corrections 
for GOODS LIRGs could be as small as the values based 
on redshifted low-z GOODS galaxies or as large as those 
based on redshifted local LIRGs. We therefore show 
these two cases as lower and upper limits. Note that 
GOODS LIRGs with L m > 10 11 L© follow the same re- 
lation, although they are slightly more symmetric than 
LIRGs with Lin > 2.5xlO n Lq. As the asymmetry 
deficit at a given redshift shows a large variation for dif- 
ferent low-z objects (see Figure [22] and Figure [23| . the 
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whole distribution of the asymmetry deficit is used to 
assign a probability for the final corrected asymmetry of 
a galaxy at a given redshift. For field galaxies, such an 
asymmetry deficit distribution is further constructed as a 
function of the galaxy i?-band brightness (see Figured]). 

4.2. Evolution of the Merger Fraction in LIRGs 

We now correct the observed asymmetry of GOODS 
LIRGs using asymmetry deficits determined from red- 
shifted low-z GOODS galaxies and local LIRGs, which 



give conservative lower- and upper-limit estimates, re- 
spectively, as discussed above. With these corrections, 
the merger fraction of GOODS LIRGs as a function of 
redshift can be derived. Note that the asymmetry cri- 
terion for mergers is usually defined as A > 0.35 (e.g. 
IConselice et al.ll2003| ). Given a systematic shift of 0.05 
between our asymmetry method and that in the liter- 
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ature (see § EJ), we adopted A > 0.3 to be consistent 6 . 
Figure O shows our result for LIRGs at Lir, > 2.5xlO n 
Lq and LIRGs at Lir > 10 11 L Q . The shaded areas are 
enclosed by lower- and upper-limits to the merger frac- 
tions. The figure shows that LIRGs are dominated by 
mergers up to z = 1.2, with merging fractions ~ 50% in 
all redshift bins. Consistent with what has been found 
in the above section, high-redshift LIRGs show slightly 
lower merger fractions than local LIRGs. 

As a test of the conclusions from the asymmetry cal- 
culations, we also visua lly classified the IR AS Revised 
Bright Galaxy Sample (jSanders et al.l l2003h as a func- 
tion of the IR luminosity. For Lir < 10 115 L Q , we used 
the Digital Sky Survey for galaxies at distance < 60 Mpc 

6 A small fraction of galaxies with A > 0.3 are not true ma- 
jor mergers, such as the highly inclined disk galaxies and some 
star-forming galaxies with multiple non-symmetrically-distributed 
bright HII regions. Here, we adopt A > 0.3 as a practical definition 
of galaxy mergers and do not correct possible contaminators. 



where the spatial resolution is < 400 pc. At 10 115 L© 
< Lir < 10 12 Lq, HST images were used. The merger 
fractions are 12%, 41%, 80% and 95%, respectively, in 
IR luminosity bins log(L m /L Q ) of [10.5, 10.99], [11.0, 
11.49], [11.5, 12.0] and [12.0, 12.4 9], where the frac- 
tion fo r the last bin was taken from ISanders fc Mirabel 
(199 6t At £ T R < 10 12 Lb, our numbers are consistent 
with ISanders fc Mirabell (|1996l ), lying between their pure 
merger fraction and merger+close-pair fraction, as some 
galaxies in close pairs have disturbed morphologies and 
are thus identified as mergers by us. Multiplying these 
fractions with the local IR luminosity functions, we ob- 
tained merger fractions of 45% and 80% at Lir > 10 11 
Lq and Lir, > 10 115 Lq as shown in Figure [6l These 
values are nearly identical to the asymmetry-identified 
ones. 

Our result of high merger fractions f or high-redshift 
LIRGs is consistent with the UDF result (|Shi et al.ll2006h 
where no corrections are applied for the observed asym- 
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metry given the much deeper exposure in the UDF. Such 
consistency further indicates that the adopted correction 
is valid for the high-redshift GOODS LIRGs. 

We discuss briefly possible cause s for the low merger 
fractions found in other studies: iBridge et al.l (|2007l ) 
identified mergers base d on asymmetr y parameters; 
IZheng et all ([20041 ) and iBell et al.l (120051) visually clas- 
sified mergers; and lLotz et aL (l2008af) used Gini— M 2 o 
to identify mergers. IBridge et alj (|2007f ) quantified the 
merger fraction of LIRGs in the Spitzer First Look 
Survey (FLS), where the imaging is shallower than 
GOODS. They did as ymmetry corrections according to 
IConselice etaU (|2005f ). which' we have found underesti- 
mate the correction for LIRGs (see § I8.3[) . This under- 
correction probably acc ounts for the l ow m erger values 
(-10-15%) they derive. IZheng et all (|2004f ) obtained a 
visually-classified merger fraction of 16% for LIRGs at 
0.4<z<1.2 in the CFRS field. Figure [20] warns of the 
bias of visual classification, as merging galaxies can look 
like normal galaxies with no or some weak asymmetric 
structures due to SB dimming. Note that Figure [201 is 
constructed for the GOODS field. At the depth of the 
images in the CFRS field, more asymmetric structures 
should be lost and high-redshift LIRGs will be artificially 
more symmetric. A simil ar bias possibly ex ists for the vi- 
sual classification used in lBell et al.l (|2005l ). The discrep- 
ancy between Gini— M20 and asymmetry may be caused 
in part by the different timcscales over which galaxies 
can be identified as mergers. For gas rich galaxy merg- 
ers, the timescale for as ymmetry is severa l times longer 
than that for Gini-M 20 (|Lotz et al.ll2008bD . 

However, we have compared the performance of CA 
and Gini-M2o on a local sample of galaxies and confirmed 
that to first order they give similar results for merger 
fractions. Since Gini-M2o is less affected by the ratio of 
S /N, it is not clear what the implications of our analysis 
would be for morphology studies using it at z ~1. Al- 



discuss these issues, a more 



would be tor morphology s 
though lLotz et all (|2008aD 
detailed investigation would be desirable. 

4.3. Evolution of the Merger Fraction in Field Galaxies 

Figure [7] shows the merger fraction for galaxies at 
M B < -19.75 and M B < -18. 94-1. 3z, compared to other 
works, where the asymmetry correction is based on the 
asymmetry deficits found for redshifted low-z GOODS 
galaxies. Our work shows a weak redshift evolution of 
galaxy merger fractions, /merger oc (l+z) 5±0 ' 3 for Mb < 
-19.75 and /merger oc (l + z) a9±0 - 3 for M B < -18.94-1.3z. 
This result is cons istent with the morphological study of 
lLotz et all ([2008al) . A gain, however, our merger fraction 
(20 -30%) is several times higher than other works except 
for IShi et all |2006). Tm g cou i^ b e cause d by the fact 
that the visual classification and asymmetric classifica- 
tion suffers f rom strong reds hift dependence as discussed 
above. Onlv lShi"et al. (2006) used images deep enough to 
detect asymmetric features as faint as for nearby galaxies 
and thus their result based on uncorrected asymmetry is 
consistent with our result. 

We notice that our relatively-high merger fraction 
(20%) at z=0.2-0.4 may require a rapid evolution to 
that at z=0. However, the current understanding of 
the galaxy merger fraction at z <0.2 is much less con- 
strained mainly due to the small volume and use of 
shallow images with poor spatial resolution. For ex- 



ample, iDe Propris et alj (|2007f ) measured a merger frac- 
tion with asymmetry for a complete galaxy sample from 
Millennium Galaxy Catalog (MGC) and found a merger 
fraction of 2-4% depending on the definition of possi- 
ble contaminators. However, the poor spatial resolution 
(1.3 arcsec) and shallow exp osure (sky noise = 26 mag 
arcsec" - 2 ) of the MGC survey ijCross et al.ll2 004). provide 
rest-frame B-band image quality for galaxies at z=0.1 
only comparable to z=l galaxies observed in the GOODS 
survey, implying a possibly significant under-estimate of 
the merger fraction due to SB effects. The image quality 
is even worse for the SDSS and 2dFGRS images. Current 
HST data do not provide a statistically large complete 
galaxy sample at z<0.2. The fraction of mergers iden- 
tifie d as galaxy pairs is als o not we ll constrained, f rom 
1% ([De Propris et al.l 120071 ) to 5% ()Lin et al.ll2008l ) for 
the same set of data but with different methods, which 
again suggests that more thorough studies are required 
to have good constraints on galaxy mergers at z<0.2. 

4.4. Infrared and B-band Luminosity Functions of 
Galaxy Mergers 

4.4.1. Methodology for Constructing IR Luminosity 
Function 

The IR luminosity functions of gala xy mergers are 
derive d broa dly follow i ng the works of iLe Floc'h et all 
1 2005^ and IShi et al.l plOSh . The 1/V max method 
(Schmidt 1968) is applied to galaxy mergers with m z < 
25, /24/im > 60 /iJy and known redshifts. The comoving 
number density in a given luminosity bin can be written 
as: 



$(Logii R )d(LogLi R ) 



(6) 



where u is the weight to correct for incompleteness and 
Knax is the maximum volume for an object to be included 
in the sample. Vmax is given by 



V m 



Zhish dV , 
il——dz, 

dz 



(7) 



where [z\ ow , z high] is the redshift over which the object 
can be detected and fl is the survey solid angle (730 
arcmin 2 ). Note that the field edge excluded for the mor- 
phology study is a small fraction (1%) of the total area. 
While z\ ow is always fixed to the low end of a redshift 
interval, the maximum redshift, zwh is defined as: 



Zhigh — mm l z bin ' Z IR ' Z " 



(8) 



where z^ 1 ^ 1 is the high end of a redshift interval, and 
z}™ 1 * is the limiting redshift at which the observed IR 
flux reaches the limiting flux of 60 /iJy. The K-correction 
for 24 /im flux is b ased on the star-forming templates of 
iRieke et all (|2009l) . z^ ;t is the limiting redshift where 

the observed z-band magnitude reaches m* mlt = 25. 
The K-correction in the z-band is based on the ACS 
photometry and KCORR ECT code (jBlanton et alj|2003t 
iBlanton fe Roweisll2007l ). 

The incompleteness correction u> includes corrections 
for the IR detection, IR objects associated with optical 
counterparts at m 2 < 25, the redshift measurement suc- 
cess rate and the criterion of secure morphologies. The 
incompleteness of the IR de tections as a function of the 
24 ^m flux density is given in lPapovich et all ((2004) . Our 
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Fig. 7. — Redshift evolution of galaxy merger fractions with given £?-band magnitude cuts. "A": mergers identified with asymmetry; 
"GM": Gini-Af20 classified mergers; "V": visually classified mergers. 



sample is limited to f2i l im = 60 /iJy at which the incom- 
pleteness is ~50%. There is no incompleteness correction 
for IR objects associated with optical counterparts with 
m z < 25. This is because m z < 25 is deep enough to 
detect nearly all optical counterparts of the IR objects 
brighter than /24 M m = 60 /iJy, given the rough correla- 
tion between th e IR lu minosity and optical luminosity in 
iLe Floc'h et al.l (|2005f ). The incompleteness correction 
for redshift measurements is defined as the ratio of the 
number of all objects to that of objects with redshift mea- 
surements within a three dimensional magnitude-color- 
color space m z -(B — i)-(V — z). Such corrections can ac- 
count for redshift measurement success as a function of 
the galaxy brightness and color. For the whole ACS pho- 
tometry sample at m z < 25, the success rate for redshift 
measurements is 68%. The final correction is for galaxies 
with secure morphologies, i.e., galaxy radius i? ga i > 15 
pixel and asymmetry uncertainty error (A) < 0.1. Such a 
correction is defined as the ratio of the number of all ob- 
jects to that of objects with secure morphologies within a 



given luminosity bin and a given redshift bin. Note that 
these corrections are small (<1.2) except for the highest 
luminosity bin in 0.4 < z < 0.6 (a correction factor of 
1.5). This indicates that the final IR luminosity function 
of galaxy mergers is mainly based on objects with secure 
morphologies. 

4.4.2. IR Luminosity Function of Galaxy Mergers 

Figure [8] shows the IR luminosity functions of galaxy 
mergers within different redshift bins. As shown in 
the figure, our I RLFs of field galaxies match those of 
ILe Floc'h et alJ (|2005h well. The 1-a uncertainties of 
the IRLFs of galaxy mergers include the Poisson noise, 
the uncertainty in the 24/im-to-LiR conversion and cos- 
mic va riance. Since ou r field size is almost as large as 
that in ILe Floc'h et al.1 (|2005l ) and 24/xm flux density is 
used to derive the total IR luminosity in both works, 
for the uncertainty in the 24/xm-to-LiR, conversion and 
co smic variance, w e simp ly adopted the result obtained 
bv lLe Floc'h~etaLl (|200l . an average of 0.2 dex upper- 
side uncertainty and 0.1 dex lower-side uncertainty. The 
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is also shown. 

final uncertainty including the Poisson noise is plotted 
in Figure [8] and listed in Table O The solid lines are 
the fit to the IR LFs of field galaxies while the dotted 
lines are the fit to the IRLFs of galaxy mergers. The 
curve for the fit is the Schetcher function. Although the 
Schetcher function fails to fit the local field galaxy IR LF 
((Perez-Gonzalez et al J 12005( 1 . our limited available data 
points do not allow us to use a formula with more pa- 
rameters, such as a double power-law. The slope of the 
Schetcher function is fixed as the average value of the fit 
to the IR LFs in the first two redshift bins. The fitted 
result is listed in Table [3l 

Figure [8] also shows the fraction of galaxy mergers as 
a function of the IR luminosity in different redshift bins. 



At all redshifts, the merger fraction increases with the 
IR luminosity, a trend sim ilar to that for local galaxies 
((Sanders fc Mirabell[l996f l. 

4.4.3. B-band Luminosity Function of Galaxy Mergers 

The rest-frame i?-band luminosity functions of galaxy 
mergers are also derived as shown in Figure [9] The re- 
sult of the merger i?-band LF is listed in Tabled] Again, 
Schetcher functions are used to fit the data points where 
the slope is fixed as the average value of the fit to the 
LF in the first three redshift bins. The fitted result is 
listed in Table [5l For the B-band luminosity functions, 
there are enough data points to break the luminosity- 
density degeneracy. The result is shown in Figure 1101 
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For the general field B-band LF, <J>*oc(l + z) 



0.1±0.5 



and 



AMg(x(— 0.4 ± 0.3)Az, whi ch are generally c onsistent 
with the literature data fsee lFaber et al.ll2007l and ref- 
erences therein). For the merger B-band LF, $*oc(l + 



z) 



0.5±0.5 



and AM * B oc(-0.4 ± 0.3)Az. 



The breaking of the luminosity-density degeneracy al- 
lows us to evaluate the evolution of the pure merger frac- 
tion: 



galaxy /merger — ^merger 



(9) 



Given r£ ga iaxy°c(l + z)°- 1±0 - 5 and n 
in the _B-band, we have 



ocfl + z \0. 5±0. 5 

merger ~ / 



B 

merger 



oc (1 



\0.4±0.7 



(10) 



The weak redshift evolution of the pure merger fraction 
/merger indicates that the observed weak evolution in a 
given i?-band magnitude is really caused by the weak 
evolution of the galaxy merger number density relative 
to the total number density. The merger rate can be 
estimated through the merger fraction and the timescale 
of 0.2-1.1 Gyr ov er which the asym metry parameter can 
identify mergers (|Lotz et al.M2008bh . 

5. DISCUSSION 



5.1. Merger-Dominated High- Redshift LIRG 
Morphologies and High Merger Fractions In 
General Field Galaxies 

While the local LIRG s are dominated by mergers 
(|Sanders fc Mirabell fi~99ll ). the result for high-redshift 
LIRGs is controve rsial, with merge r fractions from 
as low as 10-20% llZheng et all l2004t iBell et al.l 120051: 
i Bridee et all l2007t DLotz et al.l l2008aF iMelbourne et al l 
120081) to ~5 % (IShi et al.l l2006l). Motivated by the result 
of lShi et alJ ((20060 in the UDF field, we re-evaluated the 
effect of SB dimming on the morphology measurements. 
We first revised the asymmetry measurement based on 
simulations. Our new method shows a smaller scatter 
and does not suffer from a systematic offset compared 
to the one in the literature. Wc then obtained redshift-, 
IR-luminosity- and optical-luminosity-dependent asym- 
metry corrections by measuring asymmetry deficits of 
local LIRGs and low-z GOODS galaxies redshifted to 
different higher redshifts. By applying these corrections, 
we have found that high-redshift LIRGs are generally 
highly asymmetric galaxies with implied merger fractions 
~ 50% up to 2=1.2. For the general galaxy population, 
we obtained a wea k evolution of merge r fractions, consis - 
tent with r esults oflBun dv et al J (120041) iLin et all (|2008ft , 
lLotz et all (|2008al) and lde Ravel etaU (|2008f) . but with 
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Fig. 10. — Evolution in the luminosity and density of B-band lu- 
minosity functions of galaxy mergers (diamonds) and general field 
galaxies (triangles). 

several times higher merger fractions (20-30%). 

In § 14.21 and § 14.31 we summarized possible 
reasons for the lower fractions of morphologically- 
identified mergers i n other studies (|Le Feyre et alJl2000l; 
Zheng et al.l 120041: IConselice et all 120031: ICassata et al. | 
20051 iBell et all 120051: iBridge et alj l2007t lLotz et all 



2008aHMelbourne et al.ll2008f) . 

The low fraction of kinematic galaxy pairs 
(iLe Fevre et al.l l2000t iLin et alJ l2008t Ide Ravel et all 
2008) is at least partly due to the different timescales 
probed by this method compared to morphology studies. 
Also, kinematic pairs are exclusively major mergers 
while a small fraction of highly asymmetric objects 
may be minor mergers. Our high merger fraction is 
consistent with the incidence of mergers identified 
through disturbed velocity fields. Specifically, studies 
of th e velocity fields of ~60 galaxies at 0.4 < z < 
0.75 (|Neichel et al.l 120081 ; lYang et alj l200l find a low 
fraction of rotationally-supported disk galaxies and high 
fraction (41±7%) o f galaxies with complex kinematics. 
iNeichel et al.l (|2008l) also show that the auto classifica- 
tions (CA and Gini-M2o without corrections) miss half 
of the galaxies with complex kinematics and mis-classify 
them as normal disk galaxi es (also see a sim ilar result 
from numerical simulations: lLotz et al.ll2~008bl ). 

5.2. Contribution by Galaxy Mergers to the Cosmic IR 
Energy Density 

We have derived IR LFs of galaxy mergers at differ- 
ent redshifts out to z=1.2. With these IR LFs, we now 
can quantitatively evaluate the contribution of galaxy 
mergers to the cosmic IR energy density out to z ~ 1.2. 
The open circle in the top panel of Figure [11] shows the 
cosmic IR energy density estimated by integrating gen- 
er al IR LFs obtained in this wo rk, compared to the work 
of iPerez-Gonzalez et al.l (|2005l ) (grey area). The open 
star is the evolution of the cosmic IR energy density of 
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Fig. 11. — Top Panel: The total co smic IR density (open circles ) 
compared to the result (grey area) of Perez-Gonzalez et al. (2005). 
The open star symbols are the upper- and lower-limits of the cosmic 
IR density with 1-a uncertainty contributed by galaxy mergers (see 
Figure [8}. Bottom panel: The ratio and 1 — a uncertainty of the 
merger IR density to the total IR density. 

galaxy mergers obtained by integrating the IR LFs of 
galaxy mergers where the upper and lower limits corre- 
spond to the result with asymmetry corrections based on 
rcdshifted local LIRGs and low-z GOODS field galaxies, 
respectively. The bottom panel of the figure shows the 
ratio of the merger IR energy density to the total IR 
energy density. 

The galaxy mergers-produced IR energy density shows 
dramatic evolution: Sl^ 018 " 3 oc (l+z)~ 5 and (l+z)~ 6 
for the result with low and high asymmetry cor- 
rections, respectively, compared to the total IR en- 
ergy density evolution (1+ z) 3 - (|Le Floc'h et all 120051 : 
IPerez-Gonzalez et alJ 120051 ). At z > ~1 the cosmic 
IR energy density is dominated by galaxy mergers. 
This result can be expected since LIRGs dominate the 
IR energy density a t z > 0.7 (|Le Floc'h et alj 120051 : 
IPerez-Gonzalez et alj 120051) and LIRGs are dominated 
by galaxy mergers, as found in this work. 

To convert the cosmic energy IR density to the cosmic 
SFR density, the unobscured star formation as traced by 
the UV energy density needs to be included. In the local 
universe, the UV energy density corresponds to about 
30% of t he cosmic SFR density and decreases to 10% at 
z = 1.2 (|Le Floc'h et al J 120051 ). The indicated correc- 
tion from the cosmic IR energy density to the total SFR 
density for non-mergers should be larger than those for 
mergers, as the merger IR energy density has contribu- 
tions from galaxies with on-average higher IR luminos- 
ity. This implies that Figure [TT] may over-estimate the 
contribution of galaxy mergers to the total cosmic SFR 
density. 

6. CONCLUSIONS 

We present detailed morphology studies of 
intermediate-redshift (0.2< z <1.2) LIRGs and general 
field galaxies in the GOODS field by measuring galaxy 
concentration and asymmetry, with the goal to constrain 
the role of galaxy mergers in the cosmic star formation 
history. Our main conclusions are as follow: 
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(1) A new asymmetry determination method is pro- 
posed to account for the importance of underlying back- 
ground structures in accurate asymmetry measurements. 
Simulations indicate that our method does not suffer 
from a systematic offset and has small intrinsic uncer- 
tainty. 

(2) The redshift dependence of the galaxy asymmetry 
due to surface-brightness dimming is a function of the 
asymmetry itself, with larger corrections for more asym- 
metric objects. This requires careful asymmetry correc- 
tions for high-redshift galaxies. 

(3) With the necessary asymmetry corrections, high- 
redshift LIRGs are generally galaxies with high asymme- 
tries and have implied merger fractions of ~50% up to 
z=1.2, although they are slightly more symmetric than 
local LIRGs. 

(4) With similar asymmetry corrections, high-redshift 
general field galaxies show a weak redshift evolution of 
merger fractions up to z=1.2 but with a relatively high 
merger fraction (20-30%). 

(5) The B-band luminosity functions of galaxy merg- 
ers show evolution in the luminosity and density similar 
to general field galaxies. By removing the luminosity- 
density degeneracy, the pure number density of galaxy 
mergers relative to the total density shows weak redshift 



evolution, i.e., (1 + z) m with m = 0.4 ± 0.7. 

(6) The IR luminosity functions of galaxy mergers are 
derived in several redshift bins: they indicate that the 
merger fraction increases with the infrared luminosity. 
The integral of these luminosity functions shows that the 
cosmic IR energy density from galaxy mergers shows a 
dramatic redshift evolution ((1 + z)~ 5 ~ 6 ) and starts to 
dominate the total cosmic IR energy density at z >~1. 

Our study has been confined to the CA method of mor- 
phology determination. The significant corrections we 
derive for this method indicate that a similar study of 
other methods (e.g., Gini-M2o) would be desirable. 
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7. APPENDIX: NEW NOISE CORRECTION 

We revisit here the problem of the sky noise asymmetry correction. In the case of high r edshift galaxies, we fin d 
that using the minimum noise asymmetry to correct the min(^4 ga i axy+no i SO ), as proposed by iConselice et ail (2000), 
overestimates the true galaxy asymmetry and results in a relatively large error. 
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Fig. 12. — This cartoon illustrates the basic idea for rinding the real noise correction for mm(A ga i axy + no i sc ). A UDF galaxy image is 
superposed on a clean GOODS background region to create an artificial GOODS galaxy. For this galaxy, the galaxy size, concentration and 
the rotation centers giving min(yig a i axy+no i se ) can be measured. The resulting rotation centers and galaxy size are then used to measure 
the asymmetry of the clean GOODS background region, which should give the real noise correction for min(A ga i axy _|_ no ; se ). 

The noise correction is critical for an accurate measurement of the true galaxy asymmetry. For example, even for the 
minimum noise asymmetry correction (i.e., ^4 n °is C = rnin(A no i S e), referred to as the minimum method in the following), 
A™ 1 ^ is a significant fraction of A ga i axy+no i so . The median val ue of m in(A no j se ) /A ra i aX y +no is e reaches almost 70% for 
the local relatively high SB galaxy sample from lFrei et al. (1996), which IConselice et al.l (|2000[ ) used as a test sample for 
the asymmetry parameter. We have used the GOODS-South field and the UDF to test ways to measure asymmetry. 
We carried out a simulation to optimize the estimate of the true A™ 1 ^ underlying min(A ga i axy+no i SC ), using noise 
asymmetries measured in fields placed randomly around the galaxy. The basic idea of the simulation is illustrated 
in Figure [121 An UDF galaxy image was superposed on a clean GOODS background region to create an artificial 
GOODS galaxy. For this galaxy, the galaxy size, concentration and the rotation centers giving min(A ga i axy+no i se ) 
were measured. The resulting rotation center and galaxy size were then used with the galaxy removed to measure the 
asymmetry of the clean GOODS background region, which should give the true A™;g e . By putting circular regions with 
the same aperture size randomly on this background region, the randomly-produced noise asymmetry distribution was 
obtained. The goal was to find the best way to use the randomly-produced noise asymmetry distribution to estimate 
the true A™ 1 ^. In this experiment, the direct measurements of the true noise correction values let us test our method. 

Here we describe the details of the simulation. First we used the GOODS South z-band catalog to select a total of 
13 clean z-band GOODS background regions with sizes of 400x400 pixels (1 pixel=0.03"). As the GOODS catalog 
only indicates if the galaxy center is outside the region or not, we inspected each region to make sure there was no 
extended emission from galaxies with centers outside the region. These regions were distributed over the whole field 
to account for any variations over the field. We then picked 25 UDF galaxies spanning a range of S/N (50-1000) 
within half radii. The main reason for using UDF galaxies is that the UDF is deep enough to detect faint asymmetric 
structures. Our test will thus not be biased by low S/N. The UDF segmentation images were used to define the galaxy 
pixels. These galaxy pixels were superposed on each of the 13 GOODS background regions. For each UDF galaxy 
and GOODS background region, 121 artificial images were created by putting the UDF galaxy at each position of a 
11x11 grid with a separation of 10 pixels over the central 100x100 pixel regions of the GOODS background image. 
For these 121 artificial images, we measured the galaxy size and true noise asymmetry correction as described in the 
above paragraph. We then used the mean size of these galaxies to produce the distribution of random noise asymmetry 
corrections by putting circular regions randomly within the GOODS background image. Based on these randomly- 
produced corrections within the average aperture, we created 1000 random corrections for each artificial galaxy by 
multiplying the corrections within the average ape rture with the aperture size (area) of the artificial galaxy relative 
to the average aperture size (jConselice et al.l l2000f ) . In summary, for one UDF galaxy and one GOODS background 
image, we have 121 true noise asymmetry corrections corresponding to 121 artificial galaxies and for each artificial 
galaxy we have a distribution of 1000 randomly-produced noise asymmetry measures. Note that the noise associated 
with the UDF galaxy itself is not an issue in this experiment, as any noise pattern associated with the UDF galaxy 
acts like the galaxy signal in terms of searching for rotation centers giving min(A ga i axy+no i SC ). 

Figure [TBI shows comparisons between the true noise asymmetry corrections (dotted line; 121 asymmetry values in 
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Fig. 13. — Comparisons between the true noise asymmetry corrections (dotted line; 121 asymmetry values in each panel) and the 
randomly-produced noise asymmetry (solid line; 121x1000 asymmetry values in each panel) for 13 clean GOODS background regions 
superposed with a high S/N (~1000) UDF galaxy. The vertical lines are the median values of two distributions. 

each panel) and the randomly-produced noise asymmetry (solid line; 121x1000 asymmetry values in each panel) for 
the 13 clean GOODS sky regions superposed with a high S/N (^1000) UDF galaxy, while Figure [Ml shows result for a 
low S/N (~70) UDF galaxy. The vertical line indicates the median value of each distribution. Figure [TBI and Figure ITU 
indicate that using the minimum noise asymmetry to correct the noise asymmetry underlying min(j4 ga i axy+no j se ) will 
overestimate the real galaxy asymmetry in many cases. The true sky asymmetry correction actually lies approximately 
between two extreme cases, the minimum value and the median of the randomly-produced noise asymmetry. The min- 
imum case is reached when the galaxy is so faint that the min(74g a i aX y+ no iso) is actually dominated by the background. 
The median correction is appropriate when the galaxy is so bright that the min(A ga i axy+no i sc ) is dominated by galaxy 
signal, as shown in Figure [T3l Because we usually apply a lower limit S/N cut, we will rarely have the former case. 
A more important implication from these two figures is that the backgrounds always have unknown structures, which 
makes accurate noise correction impossible. Even for the same UDF galaxy and the same GOODS background image, 
the scatter of the true noise asymmetry correction is significant. For a high S/N (1000) UDF galaxy as shown in 
Figure[H the scatter is about 20% (or a (A) =0.02) and it reaches 50% (or <j(A) =0.2) for a low S/N (70) UDF galaxy 
as shown in Figure fl4l The difference in the true correction between different GOODS background regions superposed 
with the same UDF galaxy is also non-negligible. The fluctuation in the background can be caused by non-uniform 
exposure, variation in the detector response, Poisson noise and probably more important for the deep survey, a large 
number of objects below the detection limit. 

Given the significant scatter in the true noise correction caused by complicated background structures and the fact 
that we can only measure the noise asymmetry in the background region around (instead of underlying) the galaxy, 
we tried to find the best way to use the randomly-produced noise corrections to estimate the true value. To do this, 
we measured the fraction of the randomly-produced corrections below the true value for each artificial galaxy. The 
median value of the distribution of this fraction is 0.15. Therefore, the value at the 15% probability low-end tail of 
the randomly-produced noise corrections gives the best estimate of the true value. Figure [TBI compares our new noise 
correction (described as the 15% probability method in the following) to the minimum method where the minimum 
correction is simply the minimum of the randomly-produced corrections. As shown in Figure [15l the minimum-noise 
corrected galaxy asymmetry overestimates the true value, with a median offset of SA=-0.05. It may not be a severe 
problem to use the minimum noise correction as long as all the morphologies are measured in the same way. However, 
we emphasize that, as discussed above and shown in Figure [TBI Figure [T4l and Figure [TBI it is impossible to recover an 
accurate noise correction due to the complicated background structures. Figure [T5l indicates that the 68% confidence 
range of ^4„°iso " ^noiL °f our method is 0.044 compared to 0.074 for the minimum method. A scatter of 0.074 is 
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Fig. 14. — The same as in Figure IT3l but for a low S/N (70) UDF galaxy. 
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Fig. 15. — The solid histogram shows the difference between the noise correction of our 15% probability method and the true correction, 
while the dotted one shows the difference between the correction of the minimum method and the true correction. Vertical lines indicate 
the 16% and 84% low-end probability tail of these two histograms. The median offsets of our method and the minimum one are 0.00 and 
-0.05, respectively. The 68% confidence range of our method and minimum one are 0.044 and 0.074, respectively. 
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Fig. 16. — The difference (o^^ e = \A - ^notse I ) between the true correction and our 15% probabiiity method as the function of the 
standard deviation (o- ra f donl ) of the randomly-produced noise corrections. The solid line shows <T ro ^ 1 =2o- raI ? dom which can be used to 

V noise / " " noise noise 

provide conservative estimate of fjjjj for the 99% objects in the simulation. 

actually quite large given the merging criteria of A > 0.35 for the minimum method (jConselice et al.ll2003f) . Our 15% 
probability method provides almost two times less scatter in the A£°2 e - AJ^^ compared to the minimum method. 

Figure [TBI shows the difference (cr™^ = |^^^ e -^noisel) between the true correction and our 15% probability method 
as a function of the standard deviation (cr^^, om ) of the randomly-produced noise corrections. As shown above, the 
noise corrections are affected by unknown background structures and thus it would be expected that there is no 
correlation between a^olse and cr r ^^° m . However, the maximum cr^ l sc for a given t™^ 0111 is roughly correlated with 
^noisc 01 "- A s shown by the solid line, 2xcr™> s ^ om should give an upper limit to a™£l e valid for 99% of the objects in the 
simulation (e.g. effectively at about 3cr significance). 

As a summary, we propose a new noise asymmetry correction. A set of 1000 randomly-produced noise corrections 
is produced by putting circular regions in the background image around the target. The value of this distribution at 
the 15% probability low-end tail is used to correct the noise asymmetry for min(A ga i axy+no i S0 ). The error (~3tr) in 
the final measured galaxy asymmetry is taken to be two times the standard deviation of these randomly-produced 
noise asymmetries. In reality, some galaxies are always present in the field of targets. To account for this problem, we 
defined the success rate for the one-thousand circular region placements as the fraction of circular regions containing 
no galaxy signal indicated by the SExtractor segmentation image. Circular regions containing any galaxy signal were 
not used. More sets of one-thousand placements were done until one thousand successful measurements were reached. 
If the success rate for one set of placements is lower than 50%, the circular region size for the following set is taken to 
be 80% of this set. Then the measured background asym metry is reseated to th at with the original size by assuming 
background asymmetry proportional to the aperture area (jConselice et al.ll200Cfl . 

8. MORPHOLOGIES OF REDSHIFTED LOCAL LIRGS AND LOW-Z GOODS FIELD GALAXIES 

8.1. Morphologies of Redshifted Local LIRGs 

The £>-band morphologies of the local LIRGs at Lir > 10 114 L© were measured following the procedure described 
in § 12.21 Due to the high resolution and the nature of complicated morphologies, about 10% of the local LIRGs have 
multiple radii at 7y(r)=0.2. For these objects, we used the maximum radius at 7y(r)=0.2 that matches the visually- 
identified galaxy size well. 

To account for the redshift-dependence of LIRG morphologies, we also measured the redshifted B-band morphologies 
of these LIRGs at redshifts of 0.5, 0.7, 0.9 and 1.1. Since we quantified the GOODS galaxy morphologies using the 
rest-frame B-band images, we put the local LIRGs into a clean GOODS y-band sky region for local LIRGs redshifted 
to z=0.5, «-band sky region for redshifted LIRGs at z=0.7 and z=0.9, and z-band sky region for redshifted LIRGs 
at z=l.l. The exposure times for GOODS V— , i— and z— band images are around 5000, 7000 and 20000 sees, 
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Fig. 17. — The redshift evolution of the £?-band magnitude of LIRGs with Ljr > 10 11,4 Lq, where the z=0 point is for local LIRGs 
and higher-redshift points are for LIRGs in GOODS. Diamonds with error bars show the median values and 68% confidence ranges of 
the S-band magnitude distribution at 2=0, 0.4<z<0.6, 0.6<z<0.8, 0.8<z<1.0 and 1.0<z<1.2. The solid line gives the linear fit ms = 
-20.53-0.60z. 

respectively. The clean GOODS sky region has a size of 24"x24", which is large enough to cover very extended tidal 
tails. The GOODS field is too crowded to have a continuous 24"x24" clean region without any intruding object. Our 
field is obtained by merging four 12"xl2" clean regions whose sky fluctuations are consistent with each other within 
3%. 

To redshift the LIRGs, the pixel size was rebinned through IDL FREBIN.pro which rebins the pixel assuming the 
flux within a pixel is constant over the pixel area. We did not apply any size evolution, since on average the galaxy 
size R p of LIRGs does not change with redshift. The DN/sec of the local LIRG _B-band images was converted to the 
GOODS ACS counts at a given band by using their PHOTOFLAM values and then decreased by the square of the 
luminosity distance and (1+z), where (1+z) is the k-correction, since PHOTOFLAM is defined as inverse sensitivity 
in units of erg cm -2 s _1 A -1 . The resulting DN/sec is then brightened by O.6O2; magnitudes, which accounts for the 
redshift dependence of the rest-frame B-band magnitude of LIRGs as shown in Figure Q~7] 

During this pixel-rebinned and flux-rescaled process, the original LIRG image noise was correspondingly scaled 
down. To mimic the galaxy morphologies measured in the real GOODS sky region, the original galaxy background 
fluctuation should not dominate over the GOODS background in measuring the redshifted local LIRGs. To quantify 
when the scaled original galaxy background does not affect the redshifted LIRG morphology measurement, we carried 
out a test to measure the asymmetry of a clean GOODS sky region superposed on a second GOODS sky region scaled 
down by a given factor. When the second sky region has a scaled-down factor of < 0.3, its effect on the asymmetry 
measurements of the first sky region is negligible. Therefore, we measured all redshifted LIRGs with rescaled original 
background noise smaller than 0.3 of the GOODS sky noise. At redshifts of 0.5, 0.7, 0.9 and 1.1, the number of 
measurable redshifted LIRGs are 61, 83, 88 and 88 out of a total of 88 local LIRGs. Therefore, the redshifted LIRGs 
are still representative of the complete local LIRG sample. 

8.2. Morphologies of Redshifted Low-z Field Galaxies 

The galaxy sample that we used is the F-band images of 596 GOODS galaxies at 0.2 < z < 0.4 with secure 
morphology measurements. For a fiducial galaxy at z=0.3, the GOODS y-band exposure corresponds to a 10cr rest- 
frame £?-band surface brightness of 25.5 mag arcsec -2 , comparable to the depth of the HST observations of local 
LIRGs (25.0 mag arcsec -2 ). The algorithm to measure morphologies of redshifted galaxies is almost the same as for 
redshifted local ULIRGs with two differences. First we did not require the rescaled original galaxy noise to be < 0.3 
times the noise in the band where galaxies are redshifted to. Instead we used the following strategy: 

/„ow = (l-ry 2 )°- 5 B„cw + /oid (11) 
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Fig. 18. — The asymmetry deficit vs. intrinsic asymmetry. The open circles and stars are for low-redshift (0.2<z<0.4) GOODS galaxies 
redshifted to 0.4<z<0.8 with B-band luminosity evolution and no evolution, respectively. Diamonds show the evolution-independent result, 
which is the asymmetry difference using GOODS images (low S/N) and UDF images (high S/N) for galaxies with 0.4<z<0.8. The error 
bars are given at the 16% and 84% probability limits. 
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Fig. 19. — The median asymmetry and its 68% confidence range for local LIRGs and redshifted local LIRGs at z=0.5, 0.7, 0.9 and 1.1. 
The exposure time at each redshift is labelled. A linear fit to the central three points with similar exposure times gives cL4/<5.z=-0.38. 

where I ncw is the final new image, I \d is the rebinned and rescaled original galaxy, B ncw is the new background image 
in the band where the galaxy is redshifted to and 77 is the ratio of rescaled and rebinned original galaxy noise to 
the fluctuation in B ncw . This method will make sure that all redshifted galaxies will be used at all redshifts while 
maintaining the noise level. The main reason that we did not use this method for local LIRGs is that we wanted to 
quantify the effect of the complicated background structures in deep fields on the asymmetry measurements. 

Second, no additional 5-band luminosity evolution is applied, although there is a luminosity evolution in the i?-band 
luminosity function (e.g. iFaber et al1l2007t ). This is because, ideally, the asymmetry correction for a distant galaxy 
should be based on the asymmetry deficit of a redshifted low-z galaxy with intrinsically the same asymmetry and 
i?-band brightness. Applying B-band luminosity evolution will artificially underestimate the effect of SB dimming. To 
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Fig. 20. — Illustrations that apparent galaxy morphologies become less asymmetric at higher redshift as the low-surface-brightness 
asymmetric structures are lost due to surface-brightness dimming. From left to right, we show the local LIRGs and redshifted local LIRGs 
at z=0.5, 0.7, 0.9 and 1.1, respectively. The image size is about 30x30 kpc 2 . 

further demonstrate this ar gument, we measu red the asymmetry deficits for redshifted galaxies with i?-band luminosity 
evolution (AMb = -l.lAz; Fabcr et al. 2007) and no evolution. This result can be compared to the real result obtained 
by comparing the asymmetries using GOODS images and UDF images for the same galaxies. The result is shown 
in Figure [18] for galaxies with 0.4<z<0.8 within which range the UDF limiting SB in the z-band is comparable to 
the GOODS v-band limiting SB at z—0.3. As indicated by the figure, the cases with evolution underestimate the 
asymmetry correction significantly. 

8.3. The Asymmetry Deficit due to SB Dimming is a Function of the Asymmetry Itself 

Figure 033 shows the median asymmetries and the 68% confidence ranges for the local LIRGs and redshifted local 
LIRGs at z=0.5, 0.7, 0.9 and 1.1. The exposure time at each redshift is labelled. As shown in the figure, the galaxy 
appears more symmetric at high redshifts as more low-SB asymmetric structures are embedded in the background 
fluctuations due to SB dimming. A linear fit to central three points with similar exposure times (~6000secs) gives 
SA/Sz=-0.38. The asymmetry of the local LIRGs is below the extrapolation of this linear fit as the local LIRGs should 
be more asymmetric at the exposure time of ~ 6000 sees. The redshifted LIRGs at z=l.l are above the extrapolation 
of this fit as a deeper exposure at z-band dete cts more faint asymmetric structu res. Although the redshift dependence 
of the galaxy asymmetry has been noticed bv lConselice et all ([2000, 20031 l2005h . t he slope of our depen dence is much 
larger than theirs. Even compared to the test sample of local irregular galaxies (jConselice et al.|[2005l ). our slope is 
~2 times larger. We believed that the reason for such a large difference is the different local galaxy samples used to 
quantify the redshift-dependence of the galaxy asymmetry. While the test galaxy sample in Conselice's works is a 
reasonable option for study of the general high-rcdshift galaxy population, it will underestimate the galaxy asymmetry 
for studies such as quantifying merger fractions. Figure [20] illustrates how tidal tails are progressively lost and galaxies 
appear more symmetric at higher redshift. 

To further illustrate the dependence of the galaxy asymmetry deficit due to SB dimming on the asymmetry itself, 
we measured galaxy asymmetries using GOODS z-band images and UDF z-band images for about ^250 UDF galaxies 
with high S/N and large size. Figure Q] clearly shows that there is a trend of a larger asymmetry deficits due to SB 
dimming for more asymmetric galaxies. 

Figure [21] shows the asymmetry deficit for redshifted GOODS galaxies at 0.2<z<0.4. The magnitude limits cor- 
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Fig. 21. — The median asymmetry of GOODS-S galaxies at 0.2 < z < 0.4 and brighter than certain S-band magnitudes redshiftcd to 
higher redshifts. 

respond to the 70% completeness of the redshift-morphology catalog at different redshifts (see Figure H|) . Figure [21] 
further strengthens the conclusion based on local LIRGs that the redshift-dependence of galaxy asymmetry is a func- 
tion of the asymmetry itself. Brighter GOODS galaxies with higher asymmetry at z=0.3 show a steeper decline, while 
even the brightest sample shows a slower decline than the local LIRGs that have higher asymmetry. 

8.4. Asymmetry Deficits of Redshifted Local LIRGs and Low-z GOODS Field Galaxies 

Figure [5U shows the probability distributions of the asymmetry deficits for redshifted local LIRGs, where the un- 
certainty is the Poisson noise. Note that the weight has been applied for this flux limited sample to obtain a volume- 
equivalent result (see § !3.ip . The bin size for the asymmetry distribution is 0.05, the mean quadratic asymmetry error. 
The median value of the distribution in a given redshift bin can be used to correct LIRGs at that redshift. However, 
as shown in the figure, there is a significant width in the distribution, implying a variation in the asymmetry deficit 
for an individual galaxy. Such variations are expected given the complicated structures of galaxy morphologies as a 
function of the SB. 

For general galaxy populations, Figure [221 shows the probability distributions (in absolute object numbers) of the 
asymmetry deficit for low-redshift (0.2<z<0.4) GOODS galaxies redshifted to higher redshifts as functions of absolute 
B-band magnitude. The bin size for the asymmetry distribution is 0.06. At the highest redshift, the distribution for 
the faintest £>-band luminosity bin is not plotted as there is no secure morphology measurement. 
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Fig. 22. — The distribution of the asymmetry deficit for local LIRGs at Ltr, > 2.5xf0 11 L© rcdshiftcd to higher rcdshifts. 
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Fig. 23. — The asymmetry deficit distribution of low-redshift (0.2<2<0.4) GOODS galaxies rcdshiftcd to higher rcdshifts. The histograms 
are given in different B-band luminosity bins as labelled in the figure. 
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TABLE 1 

The 70% completeness cut of the 
morphology- redshift catalog shown in 
Figure [3] 



Redshift i? gal A/ rC stB MrestB 

[kpc] [mag kpc -2 ] 

0.3 2.4 -16.4 -10.3 

0.5 3.3 -17.5 -10.7 

0.7 3.9 -18.6 -11.9 

0.9 4.2 -19.1 -12.2 

1.1 4.4 -19.8 -12.4 



Note. — The apparent galaxy size, mag- 
nitude and surface brightness arc labelled in 
Figurc|3] The absolute magnitude is measured 
as M = m - 5.01og(D L ) + 5 + 2.51og(l + z), 
where L>l is the luminosity distance in pc and 
2.51og(l+z) is the K-corrcction (sec § I3.2.3L 
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Infrared Luminosity Function of Merger Galaxies 
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Note. — 1 The galaxy asymmetries are corrected based on the redshifted local LIRGs. 2 The galaxy asymme- 
tries are corrected based on the redshifted low-z (z— 0.2-0.4) GOODS field galaxies. 



TABLE 3 

Parameters of Schetcher Function Fit to IR LFs of 
Galaxy Mergers 



redshift 
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Note. — 1 The galaxy asymmetries are corrected based on the 
redshifted local LIRGs. 2 The galaxy asymmetries arc corrected 
based on the redshifted low-z {z— 0.2-0.4) GOODS field galaxies. 



TABLE 4 

_B-band Luminosity Function of Merger Galaxies 
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Parameters of Schetcher Function Fit to Merger 
B-band LFs 
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